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Abstract. 

There is significant benefit to be gained by pursuing multi-messenger 
astronomy with gravitational wave and electromagnetic observations. In order 
to undertake electromagnetic follow-ups of gravitational wave signals, it will be 
necessary to accurately localize them in the sky. Since gravitational wave detectors 
are not inherently pointing instruments, localization will occur primarily through 
triangulation with a network of detectors. We investigate the expected timing 
accuracy for observed signals and the consequences for localization. In addition, 
we discuss the effect of systematic uncertainties in the waveform and calibration 
of the instruments on the localization of sources. We provide illustrative results of 
timing and localization accuracy as well as systematic effects for coalescing binary 
waveforms. 



1. Introduction 

There is a growing realization that multi-messenger astronomy will be of critical 
importance for gravitational wave astronomy. While the concept has been discussed 
for many years (e.g. [1]), only recently has a large push towards joint observations 
with electromagnetic and neutrino detectors begun (see e.g. [21 [U [5]). Initially, 
joint observations will provide additional confidence for early detections made in non- 
stationary gravitational wave data. Later, multi-messenger observations will be critical 
in extracting the maximum scientific payoff from gravitational wave observations. For 
example, the cleanest way to demonstrate that the progenitors of short 7-ray bursts 
(GRBs) are coalescing Neutron stars would be the observation of a gravitational wave 
chirp associated with a short GRB [6]. Additionally, joint observations will likely 
provide measurements of complementary parameters, thereby breaking degeneracies 
which would exist with single messenger observations. 

The road to joint observations has already been paved, with several gravitational 
wave search results being "triggered" by external observations, such as GRBs [3 [8] 
and soft gamma repeaters j9[ [10]. More recently, work has begun to ensure 
that gravitational wave observations can be followed up by other astronomical 
observatories. In the era of regular gravitational wave observations, it is likely that 
gravitational wave alerts will be followed up by large field of view optical (such as Pan 
Starrs [11]), 7-ray (Swift [12], Fermi [13]) and radio observatories (Lofar[14l [15]) as 
well as neutrino observatories (Ice Cube [16], Antares [17]). 



Triangulation of gravitational wave sources with a network of detectors 



2 



The vast majority of electromagnetic observatories are, by their very nature, 
directional. Thus, in order for gravitational wave observations to be useful to other 
astronomers, it is necessary to extract the sky location from the gravitational wave 
signal. However, gravitational wave detectors are sensitive to signals from large 
fraction of the sky and a single gravitational wave detector provides essentially no 
directional information for a short duration source. Thus, the ability to reconstruct 
the location of a transient signal is primarily due to triangulation based on the observed 
time delays of the signal at several detectors. For more than two sites, requiring a 
consistency between the observed amplitudes will also serve to restrict the allowed 
sky positions. In particular, for three detectors, using only timing information, one 
obtains two sky locations which are mirror images with respect to the plane of the 
three detectors; generically the amplitude information can be used to break this 
degeneracy. The issue of localization of gravitational wave signals with a detector 
network has been discussed previously, and several different algorithms proposed 

us na nn on m m m cm cm m\ . 

In this paper, we consider the ability of gravitational wave detectors to localize 
transient signals by considering only the timing information available at each site. 
For concreteness, we restrict attention to elliptically polarized gravitational waves, a 
choice motivated by coalescing binary waveforms. With this simple model, we obtain 
a straightforward estimate for the expected timing accuracy and use this to evaluate 
the triangulation ability of a network of detectors. From timing information alone, 
it is only possible to measure the projection of the sky location onto the detector 
baseline. Thus for two detectors, localization to a ring in the sky is possible, while for 
three detectors a reflection degeneracy in the plane of the detectors remains. 

In a gravitational wave search, there are several sources of uncertainties which will 
affect the localization ability of the search. First, when performing a matched filter 
search, there will typically be several additional parameters in the waveform model 
(such as the masses of a compact binary). These additional parameters will serve to 
degrade the localization ability. In addition, there are likely to be differences between 
the physical waveforms and the templates used. These arise due to errors in the 
waveform family due to truncation of analytic expansions or numerical inaccuracies. 
While the inaccuracies of the waveform are independent of the detector, their effect 
on the timing accuracy will depend upon detector sensitivity. Finally, there are 
uncertainties in the calibration of the detectors. These will result in the reconstructed 
gravitational wave strain h(t) differing from the actual gravitational wave signal. 
These calibration uncertainties will have a similar effect to the use of an incorrect 
waveform template. However, the calibration inaccuracies, as well as the associated 
timing errors, will be largely independent in the different detectors. 

The layout of the paper is as follows. In Section [21 we describe the restriction to 
elliptical polarization and briefly review the coalescing binary waveform. In Section 
[3] we obtain the expected timing accuracies, and in Section 3] present the localization 
ability of the network. Finally, in Section [5l we discuss the systematic uncertainties 
and their effect on timing. Throughout, we provide an illustrative example of expected 
results for binary neutron star (BNS) and binary black hole (BBH) systems. 

2. The waveform model 

In this paper, we will focus primarily on waveforms generated during binary 
coalescence. However, much of the framework introduced is applicable to a broader 
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class of waveforms. Therefore, we begin by laying out the minimal set of assumptions 
that are made on the form of the gravitational wave, before moving on to describe the 
waveform for binary neutron star and black hole coalescences in more detail. 

We restrict attention to elliptically polarized waveforms, following a definition 
proposed by Sutton and Poprocki [28] . Specifically, a waveform is said to be elliptically 
polarized if there exists a polarization frame such that the two polarizations /i+ and 
h x of the gravitational wave are related by 

h+(t) =n + a(t) cos tp{t) , h x (t) = r] x a(t) sin if (t) , (1) 

where a(t) and 4>(t) are the amplitude and phase of the waveform and ?j+,x encode 
the relative amplitudes of the two polarizations. We make the additional requirement 
that the amplitude is slowly varying (with respect to the phase), i.e. a /a <C if. Then, 
in Fourier space, the two polarizations are related by 

h x (f)=irjh + (f). (2) 

where i] 6 [— 1, 1] provides the ratio between the two amplitudes, with unity 
corresponding to circular polarization and zero to linear polarization. 

The gravitational waveform observed at a detector can be expressed as 

h(t) = F + (9,<j>,ip)h + (t;£)+F x (6,<j>,ip)h x (t;£) (3) 

Here, F + and F x are the well known detector response functions (see e.g. |2U|) 
which depend upon the sky location (9, </>) of the system relative to the detector 
and the polarization ip. We use £ to denote any additional parameters upon which the 
waveform depends. It follows straightforwardly from © and ([3]) that the gravitational 
waveform observed in a given detector can be expressed asQ 

h(t)=Aoho(t;t) + As h«(t,$ (4) 

where 

h i (f;t)=ih (f;t). (5) 

The constants Aq and A^ depend upon the location of the source relative to the 
detector and the parameter rj introduced above. Finally, we write the waveform 
explicitly in terms of amplitude and phase as 

and h« (/; £) = iA(f; e »(/*) . (6) 

2.1. Waveforms for Coalescing Binaries 

For concreteness, let us now specialize to the waveform emitted during binary 
coalescence, where we neglect the spin of the two components. Then, the orbital 
plane will not precess and the two polarizations of the waveform can be expressed as 



•(7) 



(1/2) (1 + cos 2 i) cos20 o h (t; t a , Z? ,toi,to 2 ) 
— cos l sin 2<j) hiL (t; t a , D Q , mi, ma) 

(1/2)(1 + cos 2 l) sin20 o h (t;t o , D Q , mi,m 2 ) 
+ cos i cos 2<p (t; t Q , D Ql mi, 7712) 

Here, D is the distance at which the signal is located, D a is a fiducial distance (e.g. 1 
Mpc), t Q and O are a reference time and phase for the signal (often taken as the 

\ We make use of the two "phases" of the waveform Yiq and as these arise naturally in the context 
of coalescing binaries, as we shall see in Section l2.ll 
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coalescence time and phase) , l is the inclination angle of the binary relative to the line 
of sight and mi and m 2 are the masses of the binary's components. The two phases 
ho, hs. are the waveforms, normalized for a binary at distance D , and depend upon 
the masses of the components as well as the coalescence time t . 

The amplitude and phase of the waveform have been calculated to exquisite 
accuracy through the post-Newtonian expansion. Although the post-Newtonian 
expressions formally extend to infinite frequency, the waveform is truncated at a pre- 
specified frequency, typically the innermost stable circular orbit (ISCO). At higher 
frequencies, the finite size of the objects will cause the true waveform to differ 
significantly from the post-Newtonian expression. Furthermore, in many applications, 
the restricted post-Newtonian approximation is used where only the leading order 
amplitude term is used, while the phase is evaluated to higher post-Newtonian order. 
It is only the restricted post-Newtonian waveform which can be written in the form 
©. The amplitude satisfies A(f) oc f~ 7 ^ 6 while the detailed phasing evolution will 
depend critically upon the masses of the system [3T] . This waveform is appropriate 
for low mass binaries as the merger occurs at a higher frequency than the sensitive 
band of the detectors. We will illustrate the results in the remainder of the paper with 
numbers appropriate for a BNS signal with component masses 1.4M Q and an ISCO 
frequency of 1500 Hz. 

More recently, breakthroughs in numerical relativity have allowed for the 
calculation of entire binary black hole merger waveform [32] . Work remains ongoing 
to cover the full mass and spin parameter space. However, for non-spinning binaries 
with comparable mass components, the waveform is well understood (see, e.g. 
[33] [34], l35l l36l [37]). Indeed, a phenomenological fit to these waveforms has been 
produced in [35]. Here, the inspiral waveform extends beyond the ISCO to a merger 
frequency, after which point the amplitude evolves as A(f) oc /~ 2 / 3 , and finally 
incorporates a ringdown. For 10 — 10M© waveforms, the ISCO is at 220 Hz, but the 
phenomenological waveform continues up to 800 Hz. We will show that, by including 
the merger and ringdown information, the timing and localization accuracies for these 
waveforms can be improved dramatically. 



3. Timing Accuracy 

The parameter estimation problem has been discussed in detail in many articles. Here, 
we provide a brief overview of the method in order to fix notation, (for further details, 
see e.g. [31]). We then proceed to use the framework to address the specific problem 
of timing accuracy. In later sections, we make use of the same framework to obtain 
localization estimates and address systematic uncertainties. 



3.1. Parameter Estimation and the Fisher Matrix 

In order to decide whether there is a signal present in the data, we calculate the 
likelihood ratio of a signal h parametrized by some set fi of parameter^ being present 
in the data s, relative to the null hypothesis: 

., , p(s\h(ji)) e -(s-h(n)\s-h(n))/2 

A(/x) ~ pW) ~ ' ( ) 



§ e.g. for the coalescing binary signal introduced in section [2. II fx = (mi , m2,t ,D, 9, <j>, ip, t, (j>o). 
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where the inner product is defined as 



and S(f) is the noise power spectrum of the detector. The formalism can be used 
both for the purposes of detection and parameter estimation, as described in e.g. [40] . 
For parameter estimation, we are interested in the posterior probability distribution 
of the parameters /x, given the data s. The posterior distribution for the parameters 
fi can be obtained using Bayes' theorem as: 

where p(fi) is the prior distribution on the parameters fj, and A(/i) is the likelihood 
ratio introduced above. If we are only interested in a subset £ of the parameters fi 
and not v, we simply marginalize over the "nuisance" parameters v by integrating 
over them. 

Let us now specialise to the case discussed in Section [2] where the waveform 
is parametrized by two orthogonal phases ho and hn. satisfying ([5]), with arbitrary 
amplitudes Aq and Ail as in and substitute this waveform into the expression 
ijSJl for the likelihood. The expression is simplified by noting that the two phases are 
necessarily orthogonal, 

(*o(€)M€)) = <fcf«)IM€)> and (V£)N(S)> = °- ( u ) 

Finally, by either maximizing the likelihood with respect to Aq and Ail or by 
marginalizing over them with a uniform prior]]] we obtain 

(sMO) 2 +(sm(o) 2 

2 (ho (£)| fto(C)) 

For other parameters, it is not possible to handle the marginalization of the 
likelihood so straightforwardly. Therefore, we make use of the Fisher Information 
Matrix to obtain expected parameter estimation accuracies. Briefly, assume that 
there is a signal in the data of the form 

8 = Aof*® +A i h i (£) + n. (13) 

Further, assume that the amplitude of the signal is sufficiently large that the noise 
contribution n can be neglected. Then, expand the likelihood in the neighbourhood 
of the true signal parameter in powers of d£. Setting (ho(£)\ho(£)) = ^ we °bt am 

lnA(^)~y [1-gabdCde] (14) 
where p 2 = (Aq + A\ ) and 

g a b = (d a ho\d b ho) - (h \d a h ) (h \d b h ) - (h^\d a h )(h^\d b h ) , (15) 

is a positive definite matrix. At quadratic order, the likelihood function is 
approximated as a multi-variate Gaussian around the peak d£ = 0. The Fisher matrix 
g a t then provides an estimate of the accuracy with which the parameters £ can, in 
principle, be determined. 

|| In many cases, the sources of interest are approximately uniformly distributed in volume. This 
leads to a prior on the amplitude of A~ 4 . A uniform prior is chosen here for ease of calculation. 
For observed signals, the amplitude will be large enough that the choice of prior will not have a 
substantial effect on the parameter accuracy estimates derived later. 
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3. 2. Timing 

It is rather straightforward to utilize the formalism introduced above to investigate 
the effects of a timing error. In particular, restrict attention to the case where the 
parameter space £ is the one dimensional time parameter. In this shift of the 

time t corresponds to (frequency dependent) phase shift in the waveform and has no 
effect on the amplitude of the waveform, specifically 

h (f;t o + dt) = e 2 ^ dt h (f;t o ) . (16) 

We can obtain a posterior distribution for the timing error, making use of the 
formalism introduced above. Generally, the time will not be known apriori, certainly 
not to millisecond accuracy, so it is natural to take a uniform prior p(dt) = const. 
Then, the posterior distribution for the time offset is 

p(dt\s)ocex P ^[(h \h (dt)) 2 + {h*\h (dt)) 2 ]} . (17) 

At quadratic order, this gives0 

(h \h (dt)) 2 xtl-dt 2 [(d t ho\d t h ) - (h \d t h ) 2 ] (18) 

(h^\h (dt)) 2 ^ dt 2 {h^\d t h ) 2 . (19) 

Then, since taking the derivative of the waveform with respect to t D amounts to 
multiplication by 2nif, 

(h a \h {dt)) 2 a 1 - A7T 2 dt 2 (fho\fh ) = 1 - 4ir 2 {dt) 2 p. (20) 

Here, we have introduced /" to describe the frequency moments of the signal as 

F : = 4 r^^STjt-r- (21) 



jo S(f) 
Similarly (|19p is approximated as 

(h*\h (dt)}&2iifdt, (22) 

where / is the mean frequency defined via (|21| . Thus, at quadratic order, the timing 
distribution is 

p(dt\s) oc exp [-2p 2 n 2 a 2 f dt 2 ] where a) ^J 2 -J 2 (23) 

is the effective bandwidth of the signal. In this case, the timing estimator is unbiased 
and has a width given by 

^ = ^— ■ (24) 
2/k po~ f 

This is a simple result which encapsulates the expected timing accuracy for a given 
source. It is inversely proportional to both the signal to noise ratio (SNR) p and 
effective bandwidth o~f of the source. 

1 This expression is somewhat different that what is obtained by directly expanding Eq. dl7l . 
However, by expanding {h(dt)\h(dt)) = 1 in powers of dt, it is easy to show that the two expressions 
are equivalent. 
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Detector 


/ (Hz) 


a f (Hz) 


P = 7 


ms) 
p = 10 


Timing 
p = 7 


; accuracy (ms) 
p = 10 


Initial LIGO 
Initial Virgo 
Advanced LIGO 


150 
140 
110 


100 
140 
120 


0.23 
0.16 
0.19 


0.17 
0.11 
0.13 


0.27 
0.19 
0.21 


0.18 
0.12 
0.13 



Table 1. Timing accuracy for binary neutron stars in different detectors. The 
table gives the mean frequency and effective bandwidth of the signal. Then, the 
timing accuracies based on the quadratic approximation (1241) and the exact result 
are given for SNR 7 and 10. 



3.3. Example: Binary Coalescence 

We can apply the results obtained above to the binary coalescence waveforms of section 
12.11 As well as providing a concrete example, it will allow us to investigate the accuracy 
of the quadratic approximation used to derive Eq. (|24p above. 

We begin by considering a 1.4 — 1.4M Q BNS. The values of the mean frequency /, 
effective bandwidth a / and timing accuracy <7t are given in Table [TJ Interestingly, all 
of the detectors have similar mean frequencies and bandwidths with the broad design 
noise curve of Virgo leading to the largest effective bandwidth. Note that although 
the frequency bandwidth of Advanced LIGO is not significantly greater than initial 
LIGO, a signal at the same distance would appear with approximately twelve times 
the SNR, and therefore the timing accuracy would be fifteen times better. A detailed 
study of the timing accuracy for various different Advanced LIGO configurations has 
been investigated in [25] . 




&t(ms) 6t(ms) 

Figure 1. Left: The normalized SNR ratio against time for initial LIGO. The 
SNR for the two phases of the template are plotted, as a function of time, when 
the waveform has zero phase. The solid lines are the exact results, dotted lines 
show the SNR approximated using the quadratic approximation discussed in this 
section. Right: The timing distribution for a BNS system at a given SNR 
in initial LIGO. The dotted lines are the quadratic approximation, while the 
solid curves are the exact expression. At low SNR, the quadratic approximation 
underestimates the timing uncertainty. 

In Figure [1] we investigate the accuracy of the quadratic approximation used in 
obtaining ([2"4"| . The figure shows the SNR as a function of time for initial LIGO, 
exactly calculated and using the approximation introduced above. In both cases, 
the masses of the simulated signal and template waveform agree. The quadratic 
approximation is only good to about 0.5 ms after which it significantly underestimates 
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Waveform to ISCO 




Full waveform 






/ (Hz) 


°f (Hz) 


a t (ms) 


/(Hz) 


o- f (Hz) 


crt (ms) 


o~t (ms) 


Detector 






p= 10 






p = 10 


P/ - 10 


Initial LIGO 


120 


40 


0.37 


150 


100 


0.16 


0.14 


Initial Virgo 


90 


50 


0.34 


140 


140 


0.12 


0.10 


Advanced LIGO 


75 


50 


0.34 


120 


130 


0.11 


0.10 



Table 2. Timing accuracy based on quadratic approximation for 10 — 10A/q 
black hole binary. The results are given for the waveform truncated at the ISCO 
frequency and for the full waveform based on the phenomenological waveforms of 
[38] . We give the timing accuracy for signals at SNR 10. The final column shows 
the timing accuracy for a waveform which accumulates an SNR of 10 to ISCO, if 
found using the full waveform template. 

the recovered SNR. This leads to an underestimation of the timing uncertainty, as 
shown in the right hand plot. The approximated distributions are more sharply peaked 
than the exact ones; the timing uncertainty is underestimated by about 20% at a signal 
SNR of 7, and 15% at SNR of 10. It is only at SNR of 25 or more that the quadratic 
approximation introduces negligible error. 

Finally, consider a 10— 10M@ BBH system. Table[2]gives the effective bandwidth 
and timing accuracy for both post-Newtonian waveforms truncated at ISCO and full, 
phenomenological waveforms. Interestingly, even though only 15 to 20% additional 
SNR is accumulated after ISCO, the timing accuracy improves by as much as a factor 
of three due to the additional high frequency content of the waveform. 

4. Sky Localization from Triangulation 

Given a timing uncertainty in each of a network of detectors, this can be translated 
to a localization accuracy for the network. In the previous section, we have derived a 
simple expression for the timing accuracy for an elliptical waveform as at ~ {2tt pa f)~ l , 
where p is the observed SNR and 07 is the effective bandwidth of the signal in the 
detector. However, much of what follows makes use only of the timing accuracy, 
without reference to its derivation and would be applicable to other waveform families. 

To obtain the timing results, we have assumed that the amplitudes of the two 
phases of the waveform in each detector are independent. While this is valid for a 
single detector, for more than two detectors, these amplitudes are not independent 
since the gravitational wave has only two polarizations. This can be seen in detail 
for coalescing binaries using a simple counting argument. For three detectors we 
make 9 measurements (two amplitudes and a time in each), but these are dependent 
on only 7 parameters (£>, 6*, 0, -0, L,(j) ,t ). However, it is reasonable to assume that 
any correlations between observed detector amplitudes will only serve to improve the 
accuracies derived below. 

4.1. Two site network 

A two site network will give single measurement of timing differential and will therefore 
provide only partial localization of the signal. Suppose that the source is located at 
position R on the unit sphere, and consider two detectors separated by a distance 
(expressed in light seconds) of D. Then, the difference in the time of arrival of the 
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signal between the two sites is 

(Ti - T 2 ) = D R 



(25) 



If the two detectors have timing accuracies o~\ and o~ 2 , the distribution of the observed 
times ti and t 2 is 

(ti-Ti) 2 (t 2 -T 2 f- 



p(ti,t 2 \s) oc p(ti,t 2 )exp 



2a? 



2al 



(26) 



Localization will depend only upon the time delay {t\ — t 2 ), so we re-express (|26l) in 
terms of a fiducial arrival time and the reconstructed location r. Marginalizing over 
the arrival time, with a uniform prior, gives 

^2" 



p(r|R) oc p(r) exp 



(D-(r-R)r 

2{al + al) 



(27) 



As expected, from the timing observation in two detectors, it is only possible to 
restrict the location of the source in the direction parallel to the separation D between 
the detectors. The localization ability is improved by better timing accuracy in the 
detectors, and also by an increased baseline between detectors. 

When localizing a source, we would like to provide the smallest region of the sky 
which contains the source, with a given confidence. This requires the choice of aprior 
distribution for the sky location r. In the absence of directional mformationjjjit is 
natural to choose the prior on r to be uniformly distributed on the unit sphere, giving 



a uniform prior distribution of D • r G [—1,1] 
obtain 



Area(90%) 3.3y/a'( + . 



Then, for a 90% confidence region, we 



(28) 



4tt D 

The area is independent of the location of the signal. For the LIGO detectors D = 10 
ms so that a timing accuracy of 0.25 ms in each detector limits the signal to about 
12% of the sky. For LIGO and Virgo, the light travel time is significantly larger at 27 
ms and consequently, with the same timing accuracy, the signal can be localized to 
about 4% of the sky. 



4-2. Three site network 

The three detector result can be obtained in a similar manner. As before, we re- 
express the observed detector arrival times in terms of the reconstructed sky location 
r and fiducial arrival time to- Marginalizing over the arrival time gives 

1 



p(r|R) oc p(r) exp 



-(r-R) J M(r-R) 



where the matrix M, describing the localization accuracy, is given by 
M 



Di2Df 2 D 2 3D^ 3 DsiD^ 



a. 



(29) 



(30) 



2:i 



'31 



+ In some cases, it might be reasonable to change this assumption. For example, in many cases 
it is reasonable to restrict the prior on r to be localized to nearby galaxies |41l . Alternatively, the 
detectors' directional sensitivities make is more likely that the signal came from certain sky locations. 
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Thus M has a contribution from each pair of detectors which depends upon the 
detector separation Dy and the pairwise timing uncertainty 

2 2 

4 = ^+^ + ^ ( 31 ) 

where k ^ The timing uncertainty from a given pair of detectors is dependent 
upon the timing accuracy <7fc in the third detector. Initially, this may seem surprising, 
but arises quite naturally due to the single marginalization over the fiducial arrival 
time. Finally, we note that the two detector result ((26)) can be reproduced by taking 
er 3 -> oo. 

Since the are coplanar, M will have a zero eigenvalue and hence a degenerate 
direction e z normal to this plane. Thus, the three detector network can only restrict 
the sky location projected onto the plane formed by the detectors. In addition, since 
M is independent of the sky location — it depends solely on the location and timing 
accuracies of the individual detectors — the localization ability within the plane of 
the detectors is independent of the location in that plane. We can complete the co- 
ordinate system by introducing co-ordinates e x and e y in eigen-directions of M. In 
this basis, the sky localization distribution is 



p(r|R) oc p(r) exp 




(32) 



where (X, Y) are the co-ordinates of the source, projected onto the plane of the 
detector and (x, y) describe the recovered location. 

As before, we will use this distribution to obtain confidence regions on the sky. 
The regions will depend upon the prior distribution on r. Although a uniform 
distribution on the unit sphere does not lead to a uniform distribution on the x — y 
plane, in most cases the source localization is sufficiently accurate to treat p(x,y) as 
constant over this small region; we will make this approximation. In projecting the 
result back to the sky, we obtain two mirror sky locations (z — » —z) and an additional 
factor of (cos6>) _1 , where is the angle between the line of sight and the normal to 
the plane of the detectorsQ If we assume that the mirror degeneracy can be broken 
with amplitude consistency tests, the source can be localized with probability p to an 
area: 

Area(p) « 2na x a y [— ln(l — p)] / cos 8 . (33) 

If the reflection degeneracy cannot be broken then the area is doubled, apart from close 
to the plane of the detectors in which case the error boxes from the mirror locations 
will overlap. The best case scenario occurs when the signal is directly overhead the 
plane of the detectors. The median occurs when cos = 1/2 and gives a factor of two 
increase in localization area. 

To provide a concrete example of expected localization abilities, we consider the 
LIGO- Virgo network of detectors. The light travel time between the LIGO sites is 
10 ms, while between Virgo and the two LIGO sites is around 27 ms for both. Since 
the two LIGO detectors have similar sensitivities and a very similar orientation, a 
large fraction of sources will have similar timing uncertainties in the two detectors. 

* Of course this breaks down when the source approaches the plane of the detectors, namely 6 ~ 7r/2. 
Close to the plane of the detectors, we can approximate the sky localization by considering the extreme 
case where the source is in the plane of the detectors, and specifically at x = 1, y = 0. The uncertainty 
in the y-direction will be proportional to a y but in the z-direction, we obtain <r z oc ^/ox- 
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Therefore, we simplify by taking the timing accuracy 07 of the two LIGO detectors to 
be identical, but allow the Virgo timing a v to differ. Then, the eigen-directions of the 
matrix M are roughly aligned with the line connecting the LIGO sites e x and the line 
connecting its midpoint to Virgo i y and the localization accuracies in these directions 
are given by 



..«^L and a y ^^M±mH. (34) 
x 7ms v 22ms V ; 

For a 90% confident localization, assuming reflection degeneracy can be broken 
^M-Wfi)^) (35, 

The localization error box is double the size for the median location, cos 6 =1/2, and 
will contain two equal size, disconnected pieces if the reflection degeneracy cannot be 
broken. For the worst case scenario, where the source is in the plane of the detectors, 
we obtain an error box of over 100 deg^ for a 0.25 ms timing accuracy. For reference, 
a 1 deg 2 localization for best and median sources requires timing accuracies of 0.06 
ms, 0.04 ms in all detectors. 



4-3. Example: Binary Coalescence 

Based on the results of recent searches [32], we will take an SNR of 7 in each of the 
LIGO and Virgo detectors to be the approximate amplitude where a binary coalescence 
signal would stand above the noise background. For BNS signals of this amplitude 
Table Q] gives a timing accuracy of 0.27 ms for the LIGO detectors and 0.19 ms for 
Virgo. This gives at best case localization of 20 deg 2 . A signal would require an 
SNR of around 25 and a well located source to reduce the 90% localization ellipse 
to 1 deg 2 — certainly a possibility for the louder sources in the advanced detector 
era. The localization accuracy for 10 — 10M Q BBH waveforms is comparable to BNS, 
namely 20 deg 2 for optimally located signals at single detector SNR of 7. Interestingly, 
the inclusion of the merger and ringdown portions of the signal provide an order of 
magnitude improvement in the localization accuracy. This improvement is consistent 
with what was observed in a detailed study of parameter estimation for BBH [43] . 

The sensitive band of advanced detectors is expected to begin at around 10 Hz, 
and it it will take a BNS system over a thousand seconds to evolve from 10 Hz to 
coalescence. Thus, it is interesting to consider the possibility of localizing the source 
prior to detection to allow for early pointing of electromagnetic telescopes. The SNR 
of the signal p(t) will accumulate during the coalescence. Similarly, the frequency will 
evolve. During the inspiral phase, the frequency of the orbit evolves to leading order 
as f(t) oc (t Q — i)~ 3 / 8 where t Q is taken to as time of coalescence. Thus, making use of 
(|21[) and ([23")) we can calculate the accumulated bandwidth of the signal as a function 
of time, o~f(t). Using (|24p. we can investigate how the timing accuracy a t evolves 
during the coalescence. 

Figure [2J shows the evolution of the accumulated SNR, frequency bandwidth and 
timing accuracy over the course of the binary inspiral. While a large fraction of 
SNR has accumulated 10 seconds before merger, the bandwidth, and consequently the 
timing accuracy, is largely accrued in the last second prior to merger. Unfortunately, 
this would seem to make advanced localization of BNS systems unlikely, even with 
advanced detectors. However, electromagnetic follow-up is still worthwhile as signals 
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Figure 2. The growth of SNR p and signal bandwidth cry as a function of time 
prior to merger. The timing uncertainty at oc (ptrj) decreases as the binary 
nears merger. Although a large fraction of the SNR is accumulated well before 
merger, the bandwidth increases significantly right before merger. Hence, the 
timing uncertainty is only decreased in the last second or so prior to merger. 



associated with coalescing binaries are expected to be emitted after coalescence pH] . 
or be delayed by dispersion through the interstellar medium [45j and therefore arrive 
after the gravitational wave signal. 

5. Systematic Uncertainties 

In a gravitational wave search, there are numerous systematic errors associated to 
the waveform which affect the localization accuracy. In this section, we consider 
systematic uncertainties due to waveform and calibration errors. These issues have 
been considered in the context of detectability and parameter estimation accuracy in 
[46l l47l [48] . Errors in the template waveform might arise due to the breakdown of 
analytic approximations, or numerical inaccuracies in simulating the waveform. We 
also consider the effect of calibration errors on localization and show that they can be 
handled in a very similar manner to waveform errors. However, calibration errors are 
independent at the different detectors, whereas the waveform errors are not. 

5.1. Waveform errors 

Let us generalize the analysis introduced in section [3] by allowing for an error in the 
waveforms used for filtering. Following [49], we write: 

h (f; t ) = A(f)(l + Sa) exp[27rift + + i5<f>{f)] , (36) 

where Sa and 5<f> characterize the amplitude and phase errors in the waveform 
respectively; we continue to assume that h^(f;t Q ) ~ iho(f;t ). Since 5a and 5<f> 
are unknown functions of frequency, they cannot be treated with the standard Fisher 
matrix technology. However we can expand the likelihood in powers of 5a and 5<j) to 
second order to obtain 

P 2 

In A(5a, 5(j>, 5t) « y[l - (h (2irf5t + 5^)\h (2Trf5t + 54>)) 
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+ (h (27rf8t + 5cj ) )\h ) 2 

- (h (Sa)\h (Sa)) + (h (Sa)\h ) 2 ) . (37) 

Interestingly, there are no cross terms between the amplitude and phase errors. 
Therefore, although the amplitude errors will affect the likelihood they do not have 
an effect (at leading order) on the timing. 

Given a phasing error 8(f), it is straightforward to differentiate (|3T[) to obtain the 
timing offset St which maximizes the likelihood as 








2tt 



(38) 



In general the phasing error S(f> is not known explicitly, instead bounds on the 
maximum error S(f> m a.x are usually provided. In many cases only a maximum phase 
error is provided, independent of frequency, then the timing offset is bounded by 



■ Hf) 



2 



i"!^^. d,L wr l '- ,] )- m 

and the bound is only obtained if the phase error is maximal at all frequencies, and 
changes sign at the mean frequency /. Finally, we note that the integrand in (|39p 
gives the mean absolute deviation of the frequency. Since this is always less than or 
equal to the standard deviation, the timing offset can be bounded by 



\St\ < — 



2tt 



(40) 



For most realistic phase errors, the uncertainty obtained from (|38[) will be significantly 
smaller than (|4H)l . 

The systematic error (St from Eq. (|40|) ) and statistical error (at from Eq. ([24]) ) are 
directly comparable. Both are inversely proportional to the frequency bandwidth of 
the signal. The statistical uncertainty is also inversely proportional to the amplitude 
(or SNR) of the signal, while the systematic is independent of amplitude. Thus, for a 
phasing error of 5° , the systematic offset is guaranteed to be smaller than the statistical 
fluctuations for waveforms with an SNR less than 12. In reality, the bound in (|40f 
is rather loose and therefore a 5° phasing error is unlikely to dominate the statistical 
timing uncertainty at an SNR less than 20. 

The same waveform family will be used to search the data from all instruments. 
Therefore, the waveform error S(f> will be the same at all sites. Two detectors will 
record a different timing offset only if their power spectra (and consequently the mean 
frequency and effective bandwidth) differ. This immediately argues that between 
LIGO sites with comparable sensitivities, the timing offset due to a waveform error 
will be negligible. The timing uncertainty between LIGO and Virgo will depend upon 
the details of the waveform. Given a specific waveform and phase error, one can 
evaluate (l38l) to obtain the results. 



5.2. Calibration Errors 

Errors in the calibration of the detectors will also affect the timing accuracy. We 
have denoted the detector output or strain as s(t). This is obtained by calibrating the 
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detector output v(t) using a response junction R(f). This response function is then 
used to obtain the calibrated data:[ 



S(f) = R(f)v(f) . (41) 

In a real detector, the response function is time dependent, however we assume that 
the response function is (approximately) constant over the duration of the signal. In 
practice, the response function R c (f) which is calculated will not agree identically 
with the true response of the instrument. Using the measured response function, the 
instrumental output is calibrated to obtain the data stream: 

s c (f)=R c (f)v(f) = ^jls(f). (42) 

The errors in calibration will affect signals in the data, as well as the noise and 
consequently the noise power spectrum, S(f). 

To calculate the effect of the calibration errors, we need to evaluate the likelihood 
(fT2j) obtained when using the incorrect response function. This requires the calculation 
of both the inner product between signal and template, as well as the template 
normalization, with the incorrect response function. We will use the notation (-|-} c to 
denote the existence of calibration errors in the power spectrum used in computing 
the inner product. 

We begin by calculating the effect of calibration errors on the inner product 
between signal and template: 



(s c \h) c = 4Re J df 



Rc(f) 2 



s(f) 



4Re / d f ^TTx — = < S IM ■ (43) 



S(f) 

where we have introduced an effective waveform error due to calibration as: 

hM '(W)) k,) - ,44) 

Similarly, it is straightforward to show that calibration errors in the calculation of the 
template norm can be expressed as 

(h\h) c =(h c \h c ). (45) 

Therefore, the calibration errors can be quantified in exactly the same way as the 
waveform errors discussed in the previous section where 

(l + 5a(/))exp[i^(/)] = (^) . (46) 

This is precisely the form in which calibration errors are expressed, for example in 

HQ]. 

The timing errors for a given phase accuracy are as given in the previous section. 
However, as has already been emphasized, the calibration errors are uncorrelated 
between instruments, so expressions (|38ti40|) are directly applicable. 

J In practice, the process is generally performed in the time domain to directly produce s(t). However, 
that will not affect the discussion below as the same systematic uncertainties in calibrating the data 
still arise. 
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5.3. Example: Binary Coalescence 

Let us return once more to the coalescing binary waveforms and evaluate the effect 
of waveform uncertainties on the timing. In Figure [3] we plot the integrand of (|38[) 
for both the BNS and 10 — 10M© BBH signals. The most significant contribution to 
timing errors arise at frequencies above and below the mean frequency /, but where 
the detector still has good sensitivity. Given a model of the waveform error 6<f)(f), it 
can be integrated against the curves in Figure [3] to obtain the timing offset. 




10 2 10 3 Itl 2 I0 3 

Frequency (Hz) Frequency (Hz) 



Figure 3. The contribution that different frequencies make to the timing offset 
given in equation (J38J - The left hand plot is for the BNS system, the right 
hand plot is for the full 10 — IOMq coalescence waveform. For both the LIGO 
and Virgo detectors, frequencies far from the mean, but which still contribute to 
the signal to noise ratio, contribute most significantly. The different shapes are 
directly dependent on the instrumental noise curves. The final curve shows the 
contribution to the timing differential as a function of frequency. 



To calculate a worst case scenario for a given maximum phase offset <50 max , we 
take the absolute value of the curves in Figure [3J and integrate. The resultant errors 
for BNS are: 



\St 



LIGO 



< 



0.09ms, l^vhgol < 



5° 



0.06ms . 



(47) 



These are limits on timing errors which would result in each detector from a calibration 
error of this magnitude. The errors would only be achieved in the (rather unrealistic) 
circumstance where the phase error was +5° up to the mean frequency and then —5° 
at all frequencies above this. 

The use of an incorrect waveform family will introduce a correlated timing error 
at the different sites. In particular, it will not introduce a timing systematic between 
the two LIGO sites, since the detectors have the same noise power spectrum. Figure 
[3] shows the relative timing offset between LIGO and Virgo detectors. It can again be 
integrated to obtain the maximal timing effect of 



I^LIGO — #tvirgo| — 



0.05ms 



(48) 



For the BBH waveform, the 5° errors are slightly larger at 0.09 ms for LIGO, 0.07 ms 
for Virgo and 0.05 ms between LIGO and Virgo. 

It is illustrative to compare these systematic errors to the statistical timing 
uncertainties derived in Section[3j In the worst case scenario, a calibration uncertainty 
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Single Det 


Waveform 


Calibration 


Timin 


I (ms) 


Median 90% 


SNR 


Error (°) 


Error (°) 


LIGO 


Virgo 


Localization (deg 2 ) 


7 


- 


- 


0.27 


0.21 


40 


7 


10 


10 


0.32 


0.26 


60 


7 


20 


20 


0.45 


0.38 


120 


10 






0.18 


0.13 


17 


10 


5 


5 


0.20 


0.15 


22 


10 


10 


10 


0.25 


0.20 


36 


25 






0.06 


0.05 


2 


25 


5 


5 


0.16 


0.11 


6 


25 


10 


10 


0.19 


0.13 


20 



Table 3. Timing and localization accuracy for binary neutron star systems in 
the LIGO- Virgo network for a range of single detector SNRs and calibration 
and waveform uncertainties. The intrinsic timing uncertainties are added in 
quadrature with the calibration and waveform uncertainties. The localization 
area is calculated from timing information, with the additional assumption that 
amplitude consistency can break the reflection degeneracy. For an ideally located 
source (orthogonal to the plane of the detectors) the localization area is half the 
median value. The numbers are given for the initial/enhanced network; advanced 
detectors would yield around a 20% improvement in timing and 40% in localization 
for the same SNR. 

in phase of 15° would introduce an error equivalent to the inherent timing uncertainty 
at (single detector) SNR of 7. Since the waveform errors will not affect the timing 
between LIGO sites, they can be effectively absorbed into the Virgo timing accuracy, 
in which case 20° waveform uncertainties will produce comparable uncertainties to 
inherent timing errors at SNR 7. 

Since the waveform, calibration and statistical uncertainties in timing are 
independent, it is natural to add them in quadrature. For a signal at SNR 7 and 10° 
waveform and calibration errors, we obtain o~i < 0.32ms for LIGO and o~ v < 0.25ms 
for Virgo. This leads to a 50% degradation of the localization accuracy, so that for 
ideally located sources the area of the 90% confidence localization ellipse increases 
to 30 deg 2 , and for an average source to 60 deg 2 . Table [3] provides the localization 
accuracies for sample of SNRs and waveform/calibration uncertainties. 

6. Discussion 

We have introduced a simple method to compute the timing accuracy in a gravitational 
wave detector and, using this, derived an expression for the localization ability for 
a network of detectors. For an elliptically polarized waveform, we obtain a timing 
accuracy of a t ~ 1 / '(27r pa f) . Thus, the timing accuracy scales inversely with both the 
amplitude and "effective bandwidth" of the signal. For reference, at a signal to noise 
ratio of 7, the timing accuracy for low mass coalescing binary signals is around 0.25 
ms. This holds for both BNS and low mass BBH, although for the latter it is only by 
considering the full coalescence waveform that this accuracy can be achieved. 

For a given timing accuracy in a network of detectors, we have calculated the 
accuracy with which the source can be localized on the sky. The localization accuracy 
depends upon the timing accuracy in each of the detectors, the network geometry and 
the angle between the plane of the detectors and the signal location. A detector 
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network with widely separated detectors affords the best localization ability, and 
signals which are normal to the plane of a three detector network are localized with the 
greatest accuracy. The 90% localization ellipse for an optimally oriented source, with 
timing accuracy 0.25 ms, in the LIGO-Virgo network has an area of20deg 2 . This is 
doubled for the median location, and doubled again if the reflection degeneracy in the 
plane of the detectors cannot be broken by other considerations, such as amplitude 
consistency. 

The expressions for localization accuracy derived here could easily be extended to 
a network of more that three detectors. Although the precise form of the localization 
distribution has not been calculated, it would be similar in form to the three detector 
expression ([29)) . In concurrence with other works (e.g. [51]) this argues that additional 
detectors would most improve localization efforts by providing a large baseline between 
the new detector and existing LIGO-Virgo network, as well as breaking the reflection 
symmetry by lying well away from the plane formed by the LIGO-Virgo network. 

We have given detailed expressions for the effect of systematic errors on 
localization. The effects of waveform and calibration uncertainties is almost identical 
in a single detector. However, the waveform errors will be correlated across the 
detectors in the network, while calibration errors will not. Due to the similar sensitivity 
curves of the two LIGO detectors, waveform uncertainties will have little effect on 
localization, while they will effect the LIGO-Virgo network. Calibration errors will 
be independent at the different sites and therefore might have a larger effect on 
localization. We have obtained a bound on the timing error due to calibration errors 
as St < (<50 max /27rcr/). By comparing with the inherent timing uncertainty, we see 
that calibration uncertainty will not dominate provided <5<^ max < 1/p- Note, however, 
that only for a very specific, and unlikely, form of the calibration error will the timing 
offset be anything like this large. 

In this paper, we have not considered the effect of the mismatch between waveform 
and template parameters. This will surely degrade the localization accuracy which has 
been derived. However, as for waveform error, the similarity of the LIGO detectors' 
sensitivities means that waveform errors will have produce negligible timing effect 
between them, as has been observed in [53] ■ Furthermore, in [53][5I]> it has been 
argued that the effect of parameter uncertainties can be minimized by choosing an 
appropriate reference time. In the future, we plan to investigate the effect that 
parameter mismatch has on timing accuracies, as well as exploring in greater detail 
the effect of waveform and calibration errors for multi-detector parameter estimation. 

There are numerous simplifications and approximations which are made in this 
paper. While the basic results derived here are qualitatively correct, the detailed 
expressions for sky localization ellipses will surely be modified as these assumptions 
are relaxed. For example, if the components of the binary are spinning, then the orbital 
plane will precess during the evolution, whence the waveform will not be elliptically 
polarized. However, since the precession will be slow relative to the gravitational 
wave frequency, it is reasonable to expect that the results will be similar in the case 
of spinning binaries. Likewise, the inclusion of higher waveform harmonics [55[ 156] 
will increase the effective bandwidth of the signals; however it will also require a 
generalization of the techniques described here to correctly incorporate these signals. 

Finally, we note that the localization estimates derived here are based solely on 
triangulation between sites. Thus, we have neglected significant correlations which 
must be present in the observed signal at more than two sites. For three sites, there is 
an amplitude consistency requirement that arises from the fact that the gravitational 
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waveform has only two polarizations. Amplitude consistency should, in many cases, 
serve to break the reflection degeneracy that arises by considering triangulation alone 
and, in obtaining our most optimistic localization results, we have assumed this is 
the case. It is unlikely that amplitude consistency will further restrict the location. 
However, as has been emphasized by Searle [5 7) . for a network of three or more 
detectors, there is also a phase consistency between the observed waveform at the 
sites. For an elliptically polarized waveform observed in two sites, there is always 
an orientation and polarization such that the observed phase difference is consistent 
with a given sky location. However, for three sites, this is no longer the case and 
phase consistency can be used to further restrict possible locations. Any phasing 
requirement will naturally give a timing accuracy inversely proportional to the signal's 
mean frequency /. Thus, in the case where the mean frequency is significantly larger 
than the bandwidth, one obtains higher frequency oscillations (from phasing) on top of 
the slower falloff (from timing alone). This will lead to improvements in localization, 
and would be interesting to investigate further. However, for coalescing binaries, the 
mean frequency and bandwidth are comparable which suggests this will not be a 
significant effect. 
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